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Abstract 

A unitary transformation of the iV-ion Jaynes-Cummings Hamiltonian 
is proposed. It is shown that any approximate expression of the evolution 
operator associated with the transformed Hamiltonian retains its validity 
independently from the intensity of the external driving field. In partic- 
ular, using the rotating wave approximation, one obtains a solution for 
the iV-ion Jaynes-Cummings model which improves the standard rotating 
wave approximation solution. 



1 Introduction 

The classical Jaynes-Cummings model (JCM)0, with Hamiltonian 

Hjc = fa) a + -ujg e a z + Or (aer+ + a) <r_) , 

has become, so to say, the drosophila of quantum optics, since it gives a simple 
description of the main properties of many physical systems consisting of atoms 
in interaction with electromagnetic fields For instance, it has been for 
long at center of the attention for the study of laser cooling of ions placed 
in parabolic traps The main features of some approximate solutions of 

the JCM have been extensively studied, in particular the 'collapse and revival' 
phenomenon 
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Recently, the interest for the JCM has received a novel impulse in view of its 
applications in the fastly developing research area of quantum computing. In- 
deed, in a quantum computer (QC), information is stored in a 'quantum register' 
composed of N two-level systems representing the quantum bits, or qubits [Q. 
The storage of data and all the basic operations are implemented by inducing 
controlled dynamics on the quantum register . Since a QC is a quantum me- 
chanical system, it can perform superpositions of computation operations with a 
remarkable gain of efficiency with respect to a classical computer. For instance, 
it can be shown that the problem of factoring large numbers into primes, which 
takes exponentially increasing time on a classical computer, can be solved in a 
polynomially increasing time on a QC ||. 

Now, a couple of electronic states of an ion can be regarded as a concrete re- 
alization of the concept of qubit. Starting from this simple and effective idea, 
in 1995 Cirac and Zollcr jl(| proposed a concrete model for a ion-trap com- 
puter consisting of N atomic ions trapped in a parabolic potential well. In this 
model, each ion is regarded as a two-level system oscillating between a ground 
and an excited state (\g) and |e)), so that the number of available qubits equals 
the number of trapped ions. The control of the quantum degrees of freedom is 
achieved by addressing the ions with time, frequency and intensity controlled 
laser beams (see also |ll|). One can consider two possible arrangements of the 
laser beams: each ion individually interacting with a different tightly focalized, 
possibly zero, laser field (the original Cirac-Zollcr scheme) and the whole set of 
ions interacting with a single laser pulse [ fL2| . The mutual Coulomb interactions 
allows a communication among the trapped ions. The theoretical description of 
this physical system is given by a Hamiltonian of the Jaynes-Cummings type — 
which we will call simply the W-ion Jaynes-Cummings Hamiltonian' — namely 
by a Hamiltonian which can be reduced, via a suitable approximation, to the 
Hamiltonian of the (iV-dimensional) JCM. This kind of approximation is known 



as the rotating wave approximation (RWA) 1 1 3 . The standard RWA is essen- 



tially a perturbative method since its validity depends on the smallness of a 
parameter, the Rabi frequency, which is proportional to the intensity of the 
laser field. This is a severe limit since an intense laser field implies a fast cou- 
pling between the two internal energy levels of the trapped ions — i.e. a fast 
QC — and any operation of a ion-trap QC is performed by an accurate control 
of the dynamics of the trapped ions. 

Hence, to find a solution of the dynamical problem associated with the N- 
ion Jaynes-Cummings Hamiltonian, both simple and valid for a wide range of 
the intensity of the laser field, is of great importance for the applications. In 
the present paper, we have tried to achieve this result. The basic idea is that of 
performing first suitable transformations of the Hamiltonian and only afterwards 
applying the RWA (or any other kind of approximation). Using this strategy, 
we obtain that the transformed Hamiltonian is formally similar to the initial 
one, but the new perturbative parameter is not proportional any more to the 
intensity of the laser field but a simple bounded function of it. Thus the RWA, 
or any other approximation, once applied to the transformed Hamiltonian, will 
be scarcely sensitive of the laser intensity. In this sense, we have gone beyond 
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the standard RWA. 

The paper is organized as follows. In section the Hamiltonian describing 
the system of trapped ions in interaction with a single laser beam, expressed 
in a proper set of collective coordinates, is introduced and the standard RWA 
is resumed. Then, in section suitable transformations are performed on this 
Hamiltonian in order to obtain a 'balanced Hamiltonian' with respect to the laser 
intensity. At this point, in correspondence to the resonances of the system, one 
is able to find a good solution to the dynamical problem and this is shown in 
section ||. In section the results are extended to the general case when each 
ion interacts with a different laser beam. Eventually, a discussion of the main 
results obtained is given in section ^. 



2 The basic ion trap Hamiltonian 

A system of N equal ions of mass \x in a potential trap, with strong confinement 
along the y and z axes, and weak harmonic binding of frequency v\ along the 
s-axis (the 'trap axis'), can be described — neglecting the motion of the ions 
transverse to the trap axis — by a Hamiltonian of the following type (using 
units in which ft = 1): 

N 

^2 v p a\a v + W ({a p + a* }) , 
p=i 

where {a p } are the annihilation operators associated with a suitable system 
of collective coordinates, namely the normal coordinates of the ion chain, and 
W ({a p + a p }) is the anharmonic component of the ion- ion Coulomb interac- 
tion. Precisely, there is a nonsingular real symmetric matrix [Mj P ] (see, for 
instance, ref. Hj) such that, denoted by Xj the displacement of the j-th ion 
from its equilibrium position, the following formula holds: 

N 

x 3 =J2 M j P ( a P + a t)- (!) 

p=l 

Hence, the quantity Xj P = Mj P \a p 4- aj) can be regarded as the displacement 
of the j-th ion due to the excitation of the p-th mode. The first mode of the ion 
chain, the center of mass mode associated with the operator oi, corresponds 
classically to the ion chain oscillating back and forth rigidly. For N > 2, the 
second mode, called the breathing mode (or stretch mode), corresponds classi- 
cally to the situation in which the center of mass is motionless while each ion 
oscillates with an amplitude proportional to its equilibrium distance from the 
trap center. 

Let us now fix some j £ { 1 , . . . , N} and let us suppose that the j-th ion 
in the array is addressed by a laser beam of frequency wj, in a traveling wave 
configuration. Then, if the j-th ion is regarded as a two-level system oscillating 
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between a ground \g) and an excited state |e), the ion array will be described 
by Hamiltonian of the Jaynes-Cummings type, namely 



where 



H(t) = H + W{{a p + al})+H l (t), (2) 
N 1 

H Q :=^v p h p + -uj ge a 3 z , (3) 



P =i 

H x (t) := tt K (e^* d V* 2 + e~^ Lt o\ 2? 2 ) , (4) 

with hp = a p a p the number operator and flu = p£ the Rabi frequency, £ 
denoting the intensity of the laser field. Moreover, we have set: 

V := cxp ^ ( a v + a p]j> ( 5 ) 

where 

k L cos<j) 

Vjp : " V2^T jp (6) 

— with kj, the wavevector and <j> the angle between the x-axis and Icl — plays 
the role of a Lamb-Dicke factor. Notice that if D p {a), a e C, is a displacement 
operator associated with the p-th mode, namely 

D p (a) = exp (aa P - a* a p ) , D p (a) a p D p (a)^ — a p - a, (7) 

we have: 

V = JJ exp (i^ (a t + a p )) = n • (8) 

p— i p— i 

The Hilbert space of the total system ('pointlike' ions + internal degrees of 
freedom of the j-th ion), namely H (g> C 2 , W = L 2 (R N ), is isomorphic to — and 
will be identified with — the space TKBH, since it is often useful to work with 
operator matrices. 

Now, the dynamical problem associated with the time-dependent Hamilto- 
nian H{t) can be reduced to the dynamical problem for a time-independent 
Hamiltonian. Indeed, switching to the interaction picture with reference Hamil- 
tonian \tji L to\ and setting 
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one obtains the time-independent Hamiltonian 
H = Rt (H(t) - \u L <rf) R\ 

N 



Rt := expli^-ULtai) , (9) 



^y p h p + ^5a j z + n R (aiV^ + cj j + V 2 ^ +W, (10) 



P =i 
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where S :— uj ge —ujL is the ion-laser detuning and we have set W = W ({a p + al}) ■ 

At this point, assuming that the anharmonic term W can be neglected, the stan- 
dard rotating wave approximation consists in expanding the exponentials that 
appear in T> 2 , 



v' 2 = Yl ex p { i7 ijp (4 + a p)) 



p=i 

N ( v 2 \ 

= II y 1 + *W<4 + «p) - -f (4 + fl p) 2 + • • -J » 

then passing to the interaction picture with reference Hamiltonian 

N 



^v p h p + -5 a 



j 

2"" 
p=i 

so obtaining the interaction picture Hamiltonian 

/ N 



flirt (*) = 1 - i E ^(e*** 4 + e" 1 ^ 4 a„) + . . . e-' M a j _ + h.c. , 



and, finally, retaining only that terms in Hi nt (t) which are slowly rotating and 
at most linear in the Lamb-Dicke parameters, that in applications are usually 
assumed to be small (rjj p <C 1)- Hence, in correspondence to the three types of 
resonance 



UJ L ~0, LU ge -U L + V k ~ 0, UJg e -Ld L - V k ~ 0, 



one obtains respectively the following three types of effective interaction picture 
Hamiltonian: 



(o J _+a\^, irj jk Q, R (a\(r i + - a k a J _^j , irj jk fi R (a k <r+ - a\ ail 



The first one, does not couple the internal degrees of freedom of the ions with 
the vibrational modes, hence it describes a situation which is not useful for 
the quantum computing applications. The last one is precisely the Jaynes- 
Cummings Hamiltonian. 

We remark that the standard RWA is not always a good approximation. Its 
validity can be checked by means of a perturbative expansion of the evolution 
operator, for instance by the Magnus expansion [ p^[ . In particular, it turns out 
that it does not work in the strong field regime (large Rabi frequency SIr), since 
the Lamb-Dicke factor r}j k does not depend on the intensity of the laser field 
but only on its frequency, so that the interaction picture Hamiltonian becomes 
'too large' for a perturbative treatment. 
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3 Beyond the standard rotating wave approxi- 
mation 

In order to go beyond the standard rotating wave approximation, we will show 
that is convenient to apply it after that the Hamiltonian H has undergone 
suitable transformations. Aim of these transformations is to obtain a 'good 
Hamiltonian' which is the sum of a large component having a simple 'diagonal' 
form — i.e. such that its matrix representation in the standard basis 



{|m) ® • • • <g> \n N ) <g) \g), \m) 



\n N ) (g)\e) :n 1 ,...,n N = 1,2,...} (11) 



is diagonal — and a small spin-flip component scarcely sensitive of the Rabi 
frequency flu. 

As a first step, extending the transformation used by Moya-Cessa et al. jl6| 
for a single ion to the N-ion system, namely the transformation associated with 
the unitary operator 



Tx:= ^= (~ (V + Pt) - I (p - x>t) CT i + v ^ v \ j_ 



1 

72 



(12) 



we are able to 'linearize' the exponentials that appear in the spin-flip component 
of H (see Eq.rtlC])). Indeed, using the relations 



Ti h p T{ 



, 2 ^ p 



l P ' 4 



(13) 



(14) 



(15) 



and the fact that T\ commutes with W , one obtains the transformed Hamilto- 
nian 



where 



1 



A' 



•0 = 7", | n-lY, 1 %p v A T l = Ss * + * s > 



(16) 



So 



1 



N 

^Vph p + Q R o{ - ^6 [at + cr+ ; 

P =i 

N 

p=l 



(17) 
(18) 
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and we have dropped the unessential constant term j X^^Li rfp v v 
Let us notice that the spin-flip component of 2XJ does not depend on the intensity 
of the laser field but only on its frequency (through the Lamb-Dicke parameters) . 
On one hand, this fact is positive in the strong field regime. On the other 
hand, in the weak field regime such component cannot be regarded as a small 
perturbation. In particular, in the limit for the field intensity going to zero 
(hence: flu — > 0), one observes that f) does not assume a trivial form, as it 
happens for the Hamiltonians H and H. With this problem in mind, let us 
perform two further transformations. First, the component Sj can be easily 
'diagonalized'. Indeed, let us denote by T 2 a spin rotation by the angle 9 round 
the y-axis; namely, let us set 



cos6»/2 -sin 0/2 
sin6>/2 cos 0/2 



Now, since 



N 



where 

H := y^^pfip 
P =i 

%= -7= A2 J2 (?jp v p ( a p - a D (?- + °+) ) • 




(19) 



T 2 (V + crQ T\ = - sin 9 a{ + cos 9 (at + , (20) 
T 2 o{ T\ = cos 9 <j{ + sin 9 (V + crQ , (21) 
we have that, if the angle — it/2 < 9 < tt/2 verifies the condition 

then the following transformations hold: 

N 

v„ n,. 4- 

2 



TiSioT} = ^> p n p + ifW4 + A 2 <7i, (23) 
P =i 

„ A N 

T 2 WT} = W-i 2 Y,V jP v P (a p -al)ai 
* p=i 

1 W 

+ 8 /4+A2 E ^p^p ( a p - a l) ( a - + g +) ' ( 24 ) 

^ p— 1 

Eventually, reorganizing the various terms, we have: 

H :=T 2 f)T% = H a + H t + W, (25) 
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At this point, in the limit for the laser intensity going to zero, the small spin- 
flip component Ifj of H cancels, as it should. The large component H still 
need to be simplified, since it contains a 'nondiagonal' term (with respect to the 
standard basis). 

To this aim, let us consider the unitary operator 



T 3 := \ (V({a p }) + P({a p })t) + \ (V({a p }) V({a p }Y) a{ 

V({a p }) 
V({a p }y 



(28) 



with V({a p }) :— Yl p=1 D p (a p ), where {a p } is a set of c-numbers to be fixed in 
a suitable way. Then, observe that 



T 3 H t£ = v v ("p _ ( a p a P + a p a p ) a J z + \a p \ 2 ) + ^-VL r ^/A + A 2 o{ 
P =i 1 

A N 

+ ^ ^ Vjp Vp ^ ap ~ °^ a * ~ ^ P ~ a * p ^ ' 



— i- 



Thus, if we assume that 



a„ = i 



2V4 + A 2 



Vjp j 



we get: 



A 2 



N 



4(4 + A 2 



vjp v pi 



(29) 



(30) 



P =i 



namely the 'diagonalized' expression that we wished to obtain. Moreover, we 
have: 

T 3 (a p - at) (ai + o+) T* = (a p - o+) (V 2>t 2 + 4 # 

A >to(oi£ t2 +7> ! ). (31) 



+ i 



where 



A" 



V:=l[D p 

P =i 



VITA 2 



VITA 2 



??jp- 



Notice that, since 



lim rfo, = sign(5) rjj P , (weak field limit) 

P-R— »0 

lim % p = 0, 



(strong field limit) 



(32) 

(33) 
(34) 
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rjj p behaves as a 'laser intensity balanced' Lamb-Dicke parameter. 
Hence, we have the transformed Hamiltonian 

/ A 2 N \ — 

6 " Ts [ H + 4(4 + A 2 ) E ^ "*J T 3 = + H t + W, (35) 

where again we have dropped the unessential constant term and we have set 

N 

;y„ «.„ 4- 

2 



N 1„ 
Ho :=^2v p n p + -6a 3 z , (36) 



P =i 



with 



and 



5 := x /402 J + ,5 2 , (37) 



%= i E ^ "p K - 4) & 2 + °i ^ 2 ) 

P =i 

- lE4^(-^ t2 -4^ 2 )- (38) 
P =i 

The new Hamiltonian H enjoys very nice properties. The first point is that the 
main component Hq has a very simple diagonal form. It only differs from the 
main component of H for the constant factor 



of ai instead of 



-8 = -(LU ge -UJ L ). 



Thus, we have obtained an effective ion-laser detuning 8 which will introduce a 
correction into the standard condition for resonances (obviously, the standard 
formula is recovered in the weak field limit). 

The second — but not less important — point is that the spin-flip component 
has a good behaviour in both the weak and strong field regime; indeed: 



lim H % = (39) 



and, since limA^o Vjp/A = Vjp/^i 



. N 

R P =i 
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Moreover, if the Lamb-Dicke parameters {rjj P } are assumed to be small, rjj P <C 1, 
since 

\rjj P \ < \Vj P \ => Vj P V 2 ~ r} jp , 
1 A 1 „ 1 „ 2 

VI+A^ > 4 + A2 ^ a ? ^»a^p' 

one has that, /or any intensity of the laser field, 

. N 

&l ~ a~ fa v p ( a P ~ 4) ( a - + °+) • ( 41 ) 

p=i 

4 The evolution operator 

In this section, we will exploit the transformations introduced above in order 
to give a suitable approximate expression of the evolution operator U(H;t,to) 
associated with the ion trap Hamiltonian H: 

ij t U(H;t,t )=H(t)U(H;t,t ), U(H;t ,t ) = Id. (42) 

Parenthetically, let us notice that, if 

H(t)=H +W + e(t)H t (t), (43) 

where denotes the Heaviside theta function, i.e. if H is the Hamiltonian 
describing the situation in which the interaction with the laser field is turned 
on at t — 0, then, for to < < t, we have simply: 

U{H- 1, t ) = U(H; t, 0) exp (i(H + W)t ) ■ (44) 

First, we will assume that the ions are not too tightly packed along the trap 
axis so that we can make the approximation 

H~Ho + H x , (45) 

namely we can neglect the anharmonic component W of the ion-ion potential. 
This assumption, whose accuracy can be estimated calculating the equilibrum 
positions of the ions in the trap (see Q ) , is also necessary in order to prevent 
the ion chain from adopting a zig-zag configuration p^| . Next, we can pass to 
the interaction picture with reference Hamiltonian Hq. Then, if we set 

VS:=exp(t#o*), (46) 
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the effective Hamiltonian of the system becomes: 



JO(t): 



N 



= fa v v ( e ~^ 1 %> 4 V t - 1 at at V\ 

P =i 

+ i E [e-^t 4 a p at I>\ - * a\ o'\ V t 

P =i 

iE4^(^v:A f -e rf v + A), 



where we have set 



to 



f := u p ± S = v v ± J 4tt 2 R + (cOg e - lo l ) 2 , 



and 



V t := Yl exp (ifj jp (e~^ 1 a p + e iu * * a\)) 



(47) 



(48) 



(49) 



P =i 



In correspondence to the resonance u k ~ 0, if <C 1, applying the RWA, one 
finds that 

(50) 
(51) 



JC(t) ~ ^ ry jfe J/fe (a fe cr^ - 4 crLj 



and hence: 



tf(J7; f, t ) =s exp ( ^^(i - i ) (o fc a\ - a\ at) 



Then, since 

(a k a\ ~ 4 a -) = (-!)"' (ofc 4 a + °t + 4 a fc CT - °+) 



it follows that 

00 

U(JO;t,t ) ~ E 



(-l) m /%fe^fc 



2m! V A 



L (*-*o) 



2 m 



(a fc 4) r 



+ E 



(-!)"» (fjjkVk 



- ' (2m+l)! V A 

m— v / x 



L (*-*o) 



2m+l 



(a£a fc ) n 


-(4«fe) ro 4 







Thus, within the approximation (|50|), the evolution operator U{ JO; t, to) has the 
following form: 



cos 



,. > *>( 1 ¥ t <«-*.)VsrFT) 



afc 



3in(to(i-fo)V^) 



cos 



(%^(*-to)V^) 



(52) 
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At this point, we will use the the transformations we performed on the orig- 
inal Hamiltonian H in order to obtain the expression of the evolution operator 
U(H;t,to)- Let us resume them: 



H(t) 



Rt 



H 



Si 



H 



H ~ Hq + f/j 



where the dashed arrows denote the switching to an interaction picture. Now, 
taking care of the difference between genuine unitary transformations of the 
Hamiltonian and Dirac (interaction picture) decompositions of the evolution 
operator, we have: 

U(H;t,t Q ) = R\ exp (-i (t - t ) H^j 

= R\ T\ exp (-i (t - to) 9}) Ti 

= R\ T} T\ exp (-% (t - t ) ff) T 2 Tx 

= R\ T\ T\ T\ exp (-i (t - t ) H) T, T 2 T\ . 
Then, using the approximation H ~ Ho + we get: 



U{H;t,t c 



R\ Tl exp 



-i{t-to)Ho + H x ) T A 



R\Tlv?U(JO;t,t )T A , 



where we have set 



(53) 



with 



T3 T 2 Ti 



c a nf=i Dpfatvjp) 
a rip=i D P ( ie A vjpV n p =i D P ( ie A vjpV 



-X 



H A Y\p=i d p (KVjp + % P )/2) t n p =i d p (KVjp - Vtp)/ 2 ) 



a nf=i d p {i% P + Vw)/ 2 ) ( -4\ 
+ rr^ n tux m \/o\n^> 




± sign(A) J - 



2V4 + A 2 2 



1 



2V4 + A 2 



(55) 
(56) 



5 Extension: N laser beams 

The results obtained above can be easily generalized to the case when each ion 
in the trap interacts with a laser beam at the same time. Notice that the case 
when all the ions are interacting simultaneously with the same laser beam will 
be included as a particular case. 
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If each ion interacts with a laser beam, then the internal degrees of freedom 
of all the ions will be involved in the dynamics of the system. Hence, the Hilbert 
space of the model will be 



N 



H ® C 2W s (H © H) <8> . . . <8> (H © H), 
where H = L 2 (R N ), while the Hamiltonian will be now the following: 

N N 

2 



AT 

+ ^ (e^K^' 4 Z>? + e *«*+^) oi ©t + w, (57) 



where QP R , lo° l are respectively the Rabi and the laser frequency associated with 
the j-th laser beam (namely the one in interaction with the j-th ion), is an 
initial phase, and we have set 



v) : = n ex P (i VjP (4 + a,,)) = n ^ 0— ^ 



P =i 



P =i 



2 / 



(58) 



being a 'Lamb-Dicke matrix'. Let us perform 



on the Hamiltonian H(t) the following transformations: 

H(t) -% H -^U H = H a + H x +W~H + H x --% jC(t), 

where we recall that the dashed arrows denote a passage to an interaction picture 
and we have: 



N , . . 

Rt :=(g)ex P Q(o;it + ^)ai), 



n p =i ^(^A, %p) x Aj n p =i A>o e X ^ 



with 



1 



2j4 + A| 



± sign(Aj) 



i 

2 ±2 



(59) 

(60) 

(61) 
(62) 

(63) 
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Accordingly, the expression of the transformed Hamiltonian H — Ho + -ffj is 
given by 

iV N 

H =^2v p h p + (64) 



2 ■ 

p=i j=i 



j,p=i j 

iV 

E^.<v(<^r (65) 



j,p=i J 



with 



5j := J 40jj 2 - (w se - c^) , % p := - Aj (66) 

v ^ > yiT^f 



AT 



p=l 



^•:=n^(^)- (67) 



2 



Then, we have 



Vi :=exp(i# t) (68) 

and the expression of interaction picture Hamiltonian J7 in the general case is 
the following: 

JJ(t) = V t H : v} 

E X ^ ( e "^ P ' "/• 4 7 V ~ e^p * 4 at T>1) 
j,p=i 3 

N 

j,p=l 3 

E X %P 'V (' " ' ^ *J;t - e ^ * < A.*) , (69) 



j,P=l 

where we have set 



u t P : = v p ± h = v p ± \j + ("V - Ml) 



(70) 



and 

N 



T> j;t := J] exp (e"*"' * a p + e*"' 'a*)). (71) 
p=i 

Thus, a resonance arises in the interaction picture Hamiltonian JO if the con- 
dition u)J = is satisfied for some j,p e {!,■•■ ,N}. Notice that up to N 
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resonances may arise simultaneously, with up to N vibrational modes being in- 
volved. The evolution operator, in the RWA, is given by a tensor product of 
identity operators and single mode evolution operators like ( |52| ) (one for each 
resonant mode). 

6 Discussion 

A very frequent problem in quantum mechanics is to find an approximate ex- 
pression of the evolution operator associated with a Hamiltonian of the following 
type: 

H(t)=H + H l (t;{e k }), 

where Hq is a explicitly diagonalizablc, in the sense that it admits a complete 
set of known eigenvectors, and {£k} is a set of adjustable parameters. Typically, 
Hi(t; {£fe}) = 0, for t < 0, and the initial (t = 0) state vector of the system is an 
eigenvector of H (namely, for t < 0, a stationary state of H). Then, a pertur- 
bative treatment of the problem is possible only if Hf(t; {£k}) is 'comparatively 
small' with respect to Hq- Usually, such condition holds only for a restricted 
range of the adjustable parameters {sk}- This turns out to be a severe limit in 
many remarkable applications. 

A problem of this sort affects the standard RWA treatment of a Hamiltonian of 
the Jaynes-Cummings type describing a set of linearly trapped ions in interac- 
tion with laser beams, if the laser intensity is not small. In the present paper, 
we have shown that such problem can be bypassed. This is accomplished by 
performing suitable transformations of the Hamiltonian. These transformations 
allow us to obtain an interaction picture Hamiltonian which is stable with re- 
spect to the intensity of the laser field. Indeed, recalling formula ( |47| ) (or (|69|)), 
one observes that the dependence on the Rabi frequency is contained in the 
factors 

with < \^f]j P \ < \ \r)j p \, and 

J_ ^ 2 _ j 2 

with < |^| fjj p < j Tjj p . Hence the validity of any approximate expression of 
the evolution operator, in particular the one obtained using the RWA, will be 
scarcely influenced by the intensity of the laser field. 

There are three further points that should be highlighted. First, a simple field- 
intensity corrected condition for resonances arises in a natural way from our 
procedure, namely (in the single laser beam case): 



v p — 8 = Q, 6= JAVL 2 R + (uj ge ~uj L f. 
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Observe that the corrected detuning S is always nonnegative, differently from the 
standard ion-laser detuning 6. Thus, both positive-;? and negative-!? resonances 
will be associated with terms of the type 

a p o+, ajcrl, 

in the interaction picture Hamiltonian JC, hence formally treated in the same 
way, differently from the standard RWA. 

Second, one can check that the action on the standard basis (|ll]) of the uni- 
tary operator Ta (or T^.y) is easily computable. This means simple explicit 
expressions for the time evolution of state vectors. A detailed analysis of the 
main features of these expressions and, in paricular, of the typical 'collapse and 
revival' phenomenon will be presented in a subsequent paper 
Third, the nonharmonic component W of the ion-ion interaction is not modi- 
fied by any of the transformations performed. To show this we have retained 
it up to formula ([45|), where it is neglected. Hence its relatively simple form 
is not corrupted and it can be taken into account as a perturbation in a more 
accurate treatment. Again, such a treatment will be presented in a subsequent 
paper {§. 

Other interesting unitary transformations have been presented, with differ- 
ent aims, in the literature [jl8j jl^] |2(j] , apart from the previously cited one jl6| . 
They do not seem to share the same properties with the one presented here. 
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